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ABSTRACT 

We develop a consistent theory for dealing with transport 
phenomena in stellar atmospheres starting with the kinetic 
equations and introducing the three cases: LTE, partial LTE 
(usually called non-LTE), and non-LTE (nonlocal distribution 
functions). We present the consistent hydrodynamical eguations 
for partial-LTE, define transport coefficients, and show a method 
to calculate them. The method is based on the numerical solution 
of kinetic equations considering Landau, Boltzmann, and Focker- 
Planck collision terms. Finally we show a set of results for the 
transport coefficients derived for a partially ionized hydrogen 
gas with radiation, considering ionization and recombination as 
well as elastic collisions. The results obtained imply major 
changes in some types of theoretical model calculations and can 
resolve some important current problems concerning energy and 
mass balance in the solar atmosphere. We show that energy 
balance in the lower solar transition region can be fully 
explained by means of radiation losses and conductive flux. 
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I. INTRODUCTION 

We do not need to docuinent the enormous importance that 
transport phenomena have in stellar atmospheres, as well as in 
most areas of physics; we only mention here that thermal 
conduction is the key to understanding the solar transition 
region and that ambipolar diffusion can be very important for the 
mass and energy transport perpendicular to the magnetic field. 

By transport phenomena we mean the macroscopic transport 
associated with matter, excluding radiation. Radiation processes 
have been thoroughly treated since the beginning of astrophysical 
research. See Athay (1972) and Mihalas (1978). 

Standard transport theories have been developed in the 
context of laboratory physics for cases where radiation is 
negligible and LTE prevails (implying small gradients of physical 
parameters as well as a small radiation flux). In that context, 
the early work by Chapman, Enskog, and Brunet (hereafter CHEB) , 
as discussed by Chapman and Cowling (1936, hereafter CC), led to 
a complete compilation of transport coefficients, definitions, 
and calculation methods for non-ionized gases; the work by 
Braginskii (1965) gives the results of applying CHEB methods to 
highly ionized gases (plasmas) and some additional definitions of 

I 

frictional coefficients very usefull in laboratory experiments. 
Hochestim (1967) has made more accurate derivations of transport 
coefficients using that method. 

The case of a partially ionized hydrogen gas has been 
treated by Devoto (1966, 1968) based on CHEB methods, but these 
results again refer to laboratory conditions (particularly high 
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pressure around 1 atm). Following the same formalism, Nowak and 
Ulmschneider (1977) have calculated the thermal conductivities 
for pressure values covering the range of interest in 
astrophysics. These last three papers consider neither radiation 
nor inelastic collisions. 

We mention also the different approach followed by Spitzer 
and Harm (1953) where not only the transport coefficients but 
also the first order distribution function was numerically 
calculated for a fully ionized gas without restriction to a 
truncated polynomial expansion. 

The results of Spitzer and Harm (hereafter SH) have been 
used extensively in several fields of physics and also show the 
limits of validity of the transport coefficients. The SH results 
let to the analytical approach of Shvarts et al . (1981) for 
calculating flux limit coefficients which have been Qualitatively 
confirmed by experiments. An analytical approach applied by 
Campbell (1984) also shows a detailed calculation of the first 
order distribution function for very simplified cases. 

Here, we stress that these approaches imply the direct 
solution of a set of kinetic equations for a fully ionized gas, 
which are expressed with Focker-Planck collision terms, instead 
of the Boltzmann term as in the CHEB method. 

An approach that fits the capabilities of modern 
"supercomputers" has been suggested by Fon tenia (1985) showing . 
how one can define a set of transport coefficients and calculate 
them by numerical methods for cases of astrophysical interest 
where LTE cannot be assumed. This approach also deals with 
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partial-LTE (p-LTE) conditions and suggests a method for cases 
where even p-LTE does not apply. A similar numerical method has 
been used by Epperlein and Haines (1986) for a fully ionized gas 
with a magnetic field, but again without considering radiation 
and inelastic collisions. Another related paper by Luciani, 

Mora, and Pellat (1985) assumes a simplified kinetic equation and 
calculates an approach to the non-localized conductive flux that 
arises when p-LTE starts to fail (this is done for a fully 
ionized gas without inelastic collisions). Calculations of non- 
local effects also have been carried out by Shoub (1983), and by 
Owocki and Canfield (1986) for the solar transition region. 

In the present paper we show a set of hydrodynamical 
equations derived frcsn the kinetic equations which applies to 
stellar atmospheres without any apriori assumptions regarding the 
distribution functions or the radiation field. Thus, for plane 
parallel cases we show one can develop the kinetic equations (and 
radiative transfer equations) in a formulation with symmetrical 
and anti symmetrical parts, that characterizes 1) the LTE regime 
as having the particle and photon distribution functions close to 
Maxwell's and Planck's formulae, 2) p-LTE when only the first 
holds, and 3) non-LTE when neither applies. Then, we apply the 
method by Fontenla (1985) to compute a consistent set of transfer 
coefficients in the p-LTE, plane-parallel case, for a partially 
ionized hydrogen gas with a radiation field. 
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II. THE HYDRODYNAMICAL EQUATIONS 


These equations can be derived from the moments of the 
kinetic equations considering only particles (i.e., excluding 
radiation), and differs from some others which include radiation 
in the moments of the kinetic equations (see for example 
Anderson, 1976), The reason for the present formulation lies in 
the fact that for most cases in stellar atmospheres the radiation 
spectrum has to be solved in detail, and its agreement with 
observations is the main goal. 

Since the inelastic collision terms do not vanish when 
taking the moments of the kinetic equations, one is left with the 
moments of the Boltzmann collision term 5 for the o particles 
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where dir^ is the impulse phase-space volume element and m^ is the 
mass of particles of species o. 

With those definitions, the statistical equilibrium 
equations result in 


an 
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where n is the number density of o particles, ^ is the 
a a 

diffusion velocity, 6 is the fluid velocity (mass center 
velocity), and is the net rate of creation of o particles per 
volume unit. 



Since mass is conserved in collisions, the usual mass 


conservation equation holds 
+ V • ( p6 ) = 0 , 

where p is the mass density. 

Taking the first momemt of the kinetic equations we find 

+ V.(?) = f + f 

where 
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and is the acceleration experinced by an a particle due to 
external (or autoconsistent ) fields. Then the definition of the 
force per volume unit experienced by the whole gas of matter 
and results the net gain of particle impulse per volume unit 
due to the inelastic collisions. The last quantity equals the 
net loss of photon impulse and can be expressed in terms of the 
collisional term for radiation 


P = — f (k I ~e)n doj dv, 

c ■* V V V 

where c is the speed of light, k , I , and e have their usual 

V V V 

meaning of absorption coefficient, intensity and emissivity of 
radiation at frequency v and with direction n, and du is the 
solid angle element. 


The tensor tt contains the pressure p, the viscous 
stress f (of null trace), and the terms due to diffusion 


ir = pf + r + Z 
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where 
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and the definition w = v - (6 + ^ ). In these equations one 

Cl 

usually drops the last term because it is quadratic in V . 

a 

For stellar atmospheres , when the photon flight time over a 
characteristic length is small, from the expression of the 
kinetic equation for photons results 
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The kinetic energy equation for the matter gas is 
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with e = E e and $ = e $ and with $ being the conductive 
a a a 

a a 

energy flux for a particles given by 
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This definition of conductive flux agrees with the one from 
CC since it does not contain the thermal energy flux due to 
diffusion 



or r ^ 4 p up to first order in V . 

OL ^ OL (X 

a 

The term 6 • f is frequently dropped for non-relativistic 
cases. Again, the condition for energy balance in collisions 
gives 


^ - S' 

a 

with E being the internal energy per a particle, 
a 

There results then 
an 

e = ~ Z E ° + V»(n ) + v»(n d)l + / ( k I - e ) du dv . 
a'^at ao a -^vv v 


By using the previous equations one can easily transform the 
thermal energy equation into the entalpy equation. 
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III. THE CALCULATION OP TRANSPORT COEFFICIENTS 

Here we give an overview of the method used for the 
calculation of transport coefficients and omit the detailed 
expressions for many of the symbols used, giving only their 
definitions . 

The basic procedure is the expansion of the kinetic 
equations for the distribution function fj^ of each kind of 
particle (A means a given species with given impulse coordinates) 
into two equations, one on the symmetrical f^® and the other on 
the anti symmetrical f^® part, with respect to the direction of 
the physical parameter gradient. 

We then write the resultant set of equations in numerical 
form by assuming a discrete partition in velocity and angle and 
replace the derivative operators by finite difference expressions 
and integral operators by sums with appropriate weight factors. 

The resultant set of algebraic equations can be solved by the 
multidimensional Newton-Raphson technigue, which is equivalent 
(Fontenla 1985) to the CHEB perturbation scheme for expressing 
the distribution function. 

With the z axis parallel to the physical parameter (p,p,T,etc.) 
gradient and assuming that fp^ does not depend on x or y (plane 
parallel geometry) , we can write the kinetic equations 
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with B = v/c, where v is the velocity and y the cosine of the 
angle 6 between the impulse and the z axis. 

Since in the present formulation we assume all variables are 
only functions of one spatial coordinate, the only consistent 
magnetic field would be a homogeneous one with the only 
nonvanishing component along that coordinate. In such a case, 
for the low field values (i.e., when the Larmour radius is much 
larger than the particle free path), the field has no influence 
on the transport coefficients and was excluded by averaging the 
kinetic eguations over angle The coefficients and 

are integral functions of the distribution functions of the other 
kinds of particles and account for the Boltzmann source term and 
sink coefficient, respectively. corresponds to the external 

(or auto-consistent) force field along the z axis and can be 
written as 
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Z .e' 
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with g being the gravitational acceleration, E* the electric 
field (divided by e) along the z axis, e the proton electric 
charge, c the speed of light, Z^ the electric charge (divided by 
e), and m^^ the mass of kind A particles. 

The remaining coefficients in equation (1) are related to 
the Landau or Focker-Planck collision terms in the kinetic 
equations (Balescu 1975) 
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Pa^ is the component of the A particle impulse, and dir^ is the 
space-phase volume element of kind B particles. The customary 
definition of is 


°ij ■ , ''ab 

2m^c 


with being the components of the A particle impulse change 

due to collision with a B particle, ^AB the relative velocity, 
and da the differential cross section. 

We shall not give all the details here, but mention only 
that the photon scattering (Thompson and Rayleigh) and the 
elastic collision terms for a heavy species A produced by a much 
lighter species B are also likely to be treated as Focker-Planck 

I t 

terms giving additional contributions to 


This is 


the case for the effect on atoms due to collisions with 
electrons, a situation very difficult to treat using Boltzmann 
collision terms requiring an extremely fine partition. 

Equation (1), with definitions (2), can be expanded in its 
swmmetrical and anti symmetrical components resulting in an 
operator formulae 

O^f^ + 0®f® = n® + D^f^ + D®f® 


and 


(3) 


O^f^ + O^f® = n® + D®f® + D^f®, 
where the operators 0 are defined by 

Tz" ^ ® ) (u|-g + 

1 3 , # S, S* . — 
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Expressions (1) - (4) are also valid for photons 

I I 

provided x n' ®n null. Equation (3) for the 

Ir t' M 

radiation field is equivalent to the usual Feautrier equations 
(Mihalas 1978). 

Equation (3) can be rewritten numerically by chosing a 
partition in z, t, v, and p space, replacing the integrals by 
sums and derivatives by finite difference quotients. The set of 


numerical equations can be solved by applying standard 
techniques . 

Equation (3) covers also stationary or quasi-sta tionary 
cases in which the time derivative in (4) can be neglected and 
the resulting set of equations can be solved by the Newton- 
Raphson iteration scheme. 

In the case of p-LTE radiation, departure from the Planck 
function (i.e., the equilibrium distribution function) is not 
negligible, and the particle populations depart from the 
Boltzmann relation, but the distribution function of particles is 
close to that of Maxwell's. In this case the convergence of CHEB 
iteration scheme rules the convergence of the Newton-Raphson 
method for solving equation (3) for particles, with the radiation 
distribution function given by the radiative transfer 
equations. 

We want to stress that the p-LTE procedure described below 
cannot be applied when the particle distribution function departs 
notably from Maxwellian at some relevant velocities. 

When the p-LTE assumption is appropriate, one develops 
equation (3) up to first order in the distribution functions, 
resulting in a set of linear equations in f^r with independent 

terms Xi? ri^, where F is the zero order distribution 

9 Z r r 

function (Maxwell's function) of physical parameters (particle 
density n^, temperature T, fluid velocity U, and external force 
potential ($z)), and Xp ^nd Hp are the Boltzmann terms due to 
interaction with photons. 
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Thus, up to first order, is a linear function of physical 
parameter gradients, the external force $ and the antisymmetr ical 
part of the radiation intensity (I^®), hereafter called 
thermodynamic forces . 

From the previous hydrodynamical equations and radiative 
transfer and statistical equilibrium equations, one obtains a 
complete set of equations for solving the full problem. The 
moment cut-off procedure is acomplished by the expansion 
conditions chosen for the distribution function. In these 
equations one can express the macroscopic fluxes as some 
transport coefficients multiplied by the thermodynamic forces. 

For this p-LTE approach the coefficients have to include 
fluxes induced by radiation, and the usual transport coefficients 
depend on the anqle-averaqed radiation field intensity (mean 
intensity). In the case where LTE holds, radiation can also be 
expanded up to first order and a reduced set of transport 
coefficients can be defined. 

One important detail not usually considered (see Mihalas 
1984, p. 421) is that if consistency is required, one has to 
consider, for the absorption and emission coefficients, the 
values corrected to first order (i.e., accounting for f®). In 
this case this results in 
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V being the frequency and the distribution function for 
photons , i .e . , 


2hv 

One can then define the optical depth as 


dx 



and the synimetrical and anti symmetrical source functions 


s _ /•2hv^'» 

V " ^2 J s 


as well as the ratio 
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vdiere S® and r^ become proportional to the thermodynamic forces. 
We then find the following set of eauations 
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It should be noted that under certain circumstances, terms 
added to the usual equations can be comparable to say 

One can easily see that the present method can be extended 
to cases where a non-thermal component of some species is present 
by treating that component in a similar way as is done for 
radiation. The only condition required is that its scale of 
variation must be large compared to the mean-free-path of the 
bulk of "thermal" particles. 

Other effects which may be considered by this method are 
collective effects such as plasma turbulence. By defining 
plasmons, soli tons, etc., one can set up expressions analogous to 
kinetic equations and the corresponding interaction terms. 

One special point concerning the method we use is the 
numerical expression adopted for derivatives of fj^ with respect 
to V and y and the considerations which have to be applied to the 
corresponding limit conditions as well as the integration weight 
to be used . 

For the f® derivatives with respect to y at intermediate 
points, we have used the standard finite difference expressions, 
but for the limit points in our grid, we have assumed 

f2(w) = ry + Sy(l-y^) 

where r and s were obtained from one limit point to the next. 

Regarding the derivatives with respect to v, we have assumed 

f®(v) = g(v) F(v) 
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where F(v) is Maxwell's function and g(v) is a slowly varying 
function of v, whose derivatives are calculated at intermediate 
points by the standard finite difference quotients and at the 
limit points by assuming 


g = r v*^ + s V 


n+m 


We have taken n = 1 for the first point and n = 3 for the 
last; m was set equal to 2 for both limits (other values also 
have been tried with no significant changes noted). The integral 
weight factors were taken using the same expansion of f^. 

The last approach for derivatives and integrals follows from 
the CHEB method and was also applied by SH, This method allows 
the use of a coarse partition in v space but it affects the 
result due to the fact that the adopted form implies certain 
boundary conditions on the function f®, both at zero and at the 
high velocity cutoff. 

It is beyond the scope of the present paper to discuss the 
generality of the solutions obtained. We mention only that, as 
CC stated, accurate results are obtained in the case where the 
bulk of low-velocity particles are very close to having a 
Maxwellian distribution, and only the very high energy particles 
can depart frcsn that distribution. 

We have tested the method thoroughly by making detailed 
comparisons with the f^/F for a fully ionized plasma from our 
method and the SH results, and for the rigid-sphere gas with the 
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transport coefficients given by CC. We find remarkable agreement 
between our computed values and those from SH for low velocities, 
and we find a substantially lower value for velocities greater 
than, for instance, twice the thermal velocity. These 
differences are understandable because SH used a constant Coulomb 
logarithm and a different cutoff, which does not have any 
influence on the resulting transport coefficients. The 
comparison with the rigid-sphere gas gives agreement to within a 
few tens of percent in the thermal conductivity, which is guite 
good for our purpose. A comparison also has been made with SH 
values for the case of a Lorentz gas where there is excellent 
agreement up to velocities of three or more times the thermal 
velocity. 


IV. THE CALCULATIONS 

One of the major goals of the present study is the 
consideration of both collision terms, the Boltzmann and Focker- 
Planck collision terms (the latter in the special form due to 
Landau, 1936, for charged particles) for neutral particles, 
charged particles, and photons, considering both elastic and 
inelastic collisions. This leads us to consider a set of 
integrodif ferential equations for the ensemble of electrons, 
protons, hydrogen atoms, and photons. In the present approach, 
we do not account for excited levels of hydrogen atcan or other 
species, and we considered photo ionization and radiative 
recombination as the only inelastic processes (reactive processes 
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since they change species). At higher temperatures 
Bremsstrahlung has to be considered, and at high densities 
collisional ionization and dielectronic recombination must be 
taken into account. 

The effect of Thompson and Rayleigh scattering is small for 
the range of pressure values considered, and we treat only a two- 
interval partition in the photon spectrum for low frequencies. 
However, we have chosen a finely divided partition of freauencies 
in the range between the head of the Lyman continuum and 
approximately 285 A (a frequency 3.2 times that at the head of 
the Lyman continuum) . 

In this paper, we have adopted a grid of 20 points in v- 
space ranging from 0.227 to 5 times the thermal velocity of the 
species; we have also made calculations changing these values and 
did not find important differences. Furthermore, the plots shown 
in Figure 8 demonstrate that the derivatives and integrals of the 
functions f^ vary smoothly between the grid points. 

For the angle 0 we have considered only three values 
covering the range from 0® to 90®, and for the <|) angle we have 
taken eight values covering the range from 0® to 360®. The steps 
of the angular partitions are not small, but were taken in order 
to obtain reasonable computing time and still represent the 
angular dependence of the first order distribution functions. 

From the results obtained one sees that the functions f^ are 
nearly linear with y. This property can be used in further work 
to simplify the calculations. 
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In order to get a reasonable set of transport coefficients 
we have chosen for the medium intensity the formula 


jS ^ 2hv^ W F = W 

V 2 _/'hv \ T V 7nv \ 

c exp(pp-)-l expfp^j^Vl 

with h and k being the Planck and Boltzmann constants, 
respectively, and W and Tr the two free parameters for specifying 
the intensity and its frequency dependence (see Fontenla and 
Rovira 1985). 

The solution of the transport equations gives 
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li 





(7) 


which may be used to relate the values of to the adopted 
thermodynamic forces 2^ and Zj^, given by 


, _ dln(W) „ _ 

W dz R dz 


In the following we will use (W , Tj^) and (Z^r Zj^) to 
describe the symmetrical and anti symmetrical parts of the 
radiation intensity (I® and I® respectively). 

For the Boltzmann terms we have considered the collision 
between a pair of incident particles, named A and B, resulting in 
a pair of emergent paricles, named C and D. The resulting 
expressions are 

Xa = / ^AB ‘^"^»'c'^*C^b‘^^B 
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( 8 ) 

where (adu d* ) is the differential cross section transformed to 
c c 

the fluid rest frame. 

The cross section for electron-atom elastic collisions were 
taken from the computations of Temkin and Lakin (1961), since 
later publications are in good agreement with these earlier 
values (Bates 1962; Khan et al. 1982) and show that the errors 
are small when considering few partial waves. 

The atom-atom collisions are far more complicated and 
calculations have to take into account more than 30 partial 
waves. However, we are not interested in the fine details of 
differential cross section which arise from the different wave 
resonances, but only in the general dependence of cross section 
with angle and energy. We have used the results of Massey (1971) 
for the combined cross section 

da = dag + I dau 

where da^ and da^ correspond to the "gerade" and "ungerade" 
interaction potentials. We notice that the total cross section 
is almost constant with respect to except for the sharp 

resonance at very small velocities. We used a formula giving 
Massey's value for the total cross section, and the shape 
of a(u) was assumed to be the one corresponding to = 0.6 km 
s“^ (E = 0.004 eV). The shape of a at higher velocities may be 
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quite different from the one given by our approach, but there 
exists no detailed publication in the range of interest. At any 
rate, since the total cross section is accurate, we do not expect 
that a change in that shape will produce major errors in our 
calculations . 

In the case of the proton-atom collisions, there is also 
another physical process, the charge transfer or capture, which 
is analogous to elastic collision, but seems to be the more 
important than elastic collision in actual cases (i.e., for small 
velocities with respect to orbital electron velocity). Again, we 
have used a formula which gives the total cross section as a 
function of fitting accurately the experimental data by Fite 
et al . (1960, 1962) and having reasonable behavior at the low 
velocity limit according to the theoretical work by Smith 
(1967). Since we do not have details on the differential cross 
sections we have assumed that the shape of a(u) is given by 
Opradolce's (1984) expression for the He'^-H charge-exchange 
collisions. 

The cross sections for collisions with photons are taken 
from Allen (1962) for Thompson and Rayleigh scattering, and for 
ionization and recombination, assuming dipolar angular 
properties . 

In equation (8) we use for f^ and fg the velocities v^, and 
Vp and the angles and 0^ which are given by the energy and 

impulse conservation equations. The values of these functions do 
not correspond to the center of an interval, and for this reason 
we have scaled those values by assuming f ^ = q F as before. This 
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procedure minimizes the error, but does not elliminate it. In 
fact, when considering the Doppler displacement for atom-photon 
collisions, errors can be important for high velocity atoms. 


V. DEFINITION OF TRANSPORT COEFFICIENTS 

Frcro our calculations we obtain the coefficients for the . 

functions f^(v,y) with respect to gradients in density of all 

particles (electrons n^, protons np, and atoms n^) , mean particle 

velocity U (fluid velocity), temperature (T) , and radiation 

parameters (W, Tj^), as well as the electric field E*. (We did 

not calculate variations with respect to gravity g since they are 

usually negligible in stellar atmospheres.) 

By using these values of f^ one can calculate the flux of 

physical quantities such as electric charge (current J), )cinetic 

energy (i.e., thermal flux q.j.); one then gets a complete set of 

transport coefficients. However, the particle densities are not 

independent. In the first place, charge neutrality ( n^ = np) 

holds in the cases of interest; otherwise enormous currents 

arise. Moreover we can assume vn = vn since the difference 

e p 

would produce strong electric fields which would have an effect 

« 

2 

larger than the straight gradient difference by a factor (L/ljj) 
(where L is the characteristic length associated with density 
gradients and Ip is the Debye length); thus, the electric field 
completely masks the other force. 

A set of macroscopic parameters can be chosen such as i, p, 
U, T, W, Tj^, E*, where i = n^/n^ is the ionization ratio and p is 



the gas pressure. In this case it is customary to define the 
ambipolar diffusion velocity as the mean velocity of atoms 
with respect to the center of gravity of protons and electrons. 

Since we assume ionization is related to the other 
macroscopic parameters through statistical eouilibrium, the 
degree of ionization is not an independent variable, and 


Z. 
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2(n, 


"a' 




- Zp + + {1 + 




where 


_ dln(i). „ _ dln(T). „ _ dln(P) 

i ~ 6z ' T dz ' p " dz 

_ dln(W) . _ - i dU 

W " dz ' R ■ dz ' U ■ c dZ 

and where n = (n^ + n_ + n_), Xt, = hv /kT„, is the frequency 

“ p » O K O 

at the head of the Lyman continuum. 

This results in a smaller set of independent transport 
coefficients that relate the fluxes (and, of course, the 
corrections in the radiative transfer equations, r and S ) to 

V V 

the thermodynamical forces Zp, Zy, Z.p, Z^, Zp, E*. 

Vie called these independent transport coefficients 
fi*p» where the asterisk has to be replaced 

by the symbol corresponding to the flux (i.e.. A, J, and T for 
the ambipolar velocity, electric current, and thermal flux, 
respectively) . 
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It is customary to define the electric conductivity 
as a = “JJjg and to define the coefficients r^ = - to 

give the ratio of the electric field E^* which will result in the 
zero electric current produced by the force . According to . 
this definition. Eg* is the value of the electric field which is 
naturally achieved when the Zg force is applied and after the 
electric charge has been redistributed to some equilibrium state, 

and rg = Eg*/Zg. 

Customary definitions regarding ambipolar diffusion are 
confusing and do not apply here because they refer to non 
reacting species or to cases where the reaction velocity is very 
small, which is not valid for hydrogen in stellar atmospheres. 

We introduce ambipolar diffusion coefficients Dg , which are 
of importance in astrophysics, by imposing only the condition of 
null electric current. Then 


D = n. + 

a Aa 


^AE ^a' 


The other transport coefficients of great importance in 
astrophysics are those related to the thermal conductive flux q, 
which are customarily defined by subtracting the entalpy flux 
frcxn the kinetic energy flux q-j.. In this case, conditions 
regarding the anullation of ambipolar velocity are not 
applicable, and the only consistent restriction is the one 
regarding electrical equilibrium. 

The definition of the thermal conductivity \ leads to the 


expression 
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AT - (2 kT ^nE ” ^TE ^ 

where n are the coefficients corresponding to the particle 

Ho 

flux . 

If more species are to be considered, one can define their 
diffusion velocity with respect to the center of mass, as 
mentioned above, and all the other coefficient definitions will 
not be altered since the times involved in reaching the zero 
diffusion concentration equilibrium are extremely large compared 
to the characteristic times for other physical processes to occur 
in stellar atmospheres. 

Note that the ambipolar diffusion velocity is not one of the 
species-diffusion velocities defined previously but forms, with 

the electric current, a complete and independent set from which 

the species-diffusion velocities can be derived 

V = - - J 

e (2 + Y ) y) 

V = - ^ 

p (2 + Y ) (1 + y) 

I 

V = V r 

a A ^2 +y' 

where 


Y 
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When replacing the above formulae in the hydrodynamical 
equations, the usual Joule heating term appears, as well as the 
internal energy transport due to ambipolar diffusion. 


VI . RESULTS 

Using the above definitions, we have calculated the 
transport coefficients for a grid of physical conditions 
characterized by the values of p, T, W and Tj^ (and assuming U = 
0). In addition we have also calculated the coefficients for Zu 
but since they are of litle interest in most cases, we do not 
show the results here. 

Since all the combinations of values for the physical 
parameters lead to a large number of possibilities, we have 
chosen only three sets of values for W and Tj^. The first set 
corresponds to optically thin matter (in the Lyman continuum) 
lying over the solar chromosphere, such as in solar prominences 
(W = 0.0022, Tj^ = 8000 K); the second set accounts for the 
enhanced radiation field above solar active regions (W = 0.022, 

Tj^ = 8000 K); the third set applies to the case of LTE (W = 1, Tj^ 
= T). 

Table 1 illustrates the grid of values for p and T used in 
our calculations. In all cases we obtained the ionization from 
the statistical equilibrium equation neglecting collisional 
ionization and recombination. We did not include stimulated 
emission, which is negligible in the cases considered here, since 
T << 157,600 K. 
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Neglecting three-body collisions in ionization and 
recombination is a good aproximation in the present cases (with 
low pressure values); we have checked the corresponding rates and 
they are negligible for p < 1000 dyne cm“^. 

The b]^ value was obtained form the expression 

E^(X) - WE^(X + Xj^) 

" W E^(Xj^) 


with X = hv^/kT, which for T and Tj^ <<157,600 K can be 
approximated by neglecting stimulated emission and by replacing 
Ejfx) by x'^ e“^.. Then 


X 

= e 

WTRe 


R 

X • 


In the following, values of transport coefficients will be 
expressed in cgs (ESU) units. We have also defined some standard 
transport coefficients as reference values. For the electrical 
conductivity we have chosen Braginskii’s (1965) expression for a 
fully ionized plasma, with InA = 10, 

Og = 1.4 X 10^ T^/^. 


For the thermal conductivity we have chosen the SH value 
given by Allen (1962) (based on InA = 10), 

Ag = 10“® T^/^. 
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In the case of the thermoelectric coefficient r-j., we have 
used the dependence on temperature given by SH, 

r^ = 10^ T . 
s 

The standard coefficients given above are provided by a 
well-established theory for fully ionized hydrogen plasmas. 
However, the remaining coefficients do not have such standard 
reference values, so we have chosen these coefficients in order 
to obtain a smooth variation at the high ionization limit. For 
the the rmo-ambi polar coefficient D.J., we estimate the asymptotic 
behavior as 


= 30 T^. 
s 

It is difficult to summarize all the results so we only show 
the data we believe to be of major importance in stellar 
atmospheric modeling. 

The ratio of the electrical conductivity to its standard value 
( ff/cr ) is shown in Figures la,b. Note that in Figures la and lb 

w * 

(solar-type and enhanced solar-type radiation fields), the ratio 
is between 0.4 and 0.8, but in the LTE case shown in Figure Ic , 
it decreases sharply for T < 6000 K (one must remember that for 
normal astrophysical abundances, because of the contributions of 
the other elements, the ratio ne/n never drops below 0.0001 which 
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is the smallest in our calculations). At moderate to hiqh 
ionization the variation of the ratio with temperature and 
pressure is due to the variation of the Coulomb logarithm 
(In A). From these results one can estimate the electric 
conductivity which has to be used when dealing with MHD theory 
for photospheric or subphotospheric layers of cold stars where 
magnetic diffusion and Joule heating can be severely 
underestimated when using Spitzer or Braginskii values. 

The thermoelectric coefficient r.p relative to its standard 
value is plotted in Figures 2. For medium to high ionization , the 
ratio is almost constant between 0.5 and 0.6. For cases of very 
small ionization in the LTE case the ratio also goes down very 
steeply, appearing to have a sharp minimum which is more 
pronounced for low pressure reaching negative values at some very 
low ionization, as shown in Figure 2c. 

The thermal conductivity (X) is shown in Figure 3. Note in 
Fig. 3a that (even at medium ionization) the thermal conductivity 
can rise to around 40 times the Spitzer value for T = 5000 K, but 
that it is only slightly higher than the usual value for 
T = 50,000 K. 

Figure 3c shows the asymptotic behavior at low ionization • 
(corresponding to the atom heat flux) and the typical behavior of 
the fully ionized plasma at the higher ionization. (Again the 
dependence with temperature and pressure is due to the variation 
of the Coulomb logarithm.) 

Figure 4 shows the shape of the guotient (D.p/D.p ) . Parts a 
and b show that the quotient is almost constant and inversely 
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proportional to the pressure, except for the case with higher 
pressure and lower radiation field. However, Fig. 4c shows 
almost constant values at high ionization that are nearly 
inversely proportional to pressure, but at low ionization and low 
pressure, the coefficient reduces and and even changes sign. 

Radiation-related coefficients, in turn, are much harder to 
describe and are not defined by previous theory. Because of the 
different usage when applying the two types of radiation fields, 
we have plotted the coefficients regarding the 2^^ force for parts 
a and b and those regarding the Zj^ force for parts c of Figures 
5, 6, and 7. The standard values were taken as 






20 





From our calculations some numerical noise arises in case of 
low pressure and high radiation field. This noise can be traced 
to the fluctuations in the lower density species functions f^(v) 
and results frcxn the limited numerical guadrature used in the 
calculations, especially for the photon spectrum and angular 
variables; usually the fluctuation increases with velocity. One 
source of this "noise," the Doppler effect on the inelastic 
collision terms, poses a difficult problem. Here, we simply do 
not consider the radiative coefficients above the limit where 
fluctuating values start. 
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Figure 5a and 5b show the shape of (r^/r^ ), and Figure 5c 
shows ( Tj^/rp ). The former exhibits an almost constant value for 
the higher pressure and lower radiation field and a behavior that 
resembles the Dp curves. Note, however, that the higher pressure 
and lower radiation intensity is the only case where the r^ 
coefficient is positive, being lower than zero in the others. 
Figure 5c shows, however, that for high ionization the ratio is 
small and inversely proportional to pressure, but is negative and 
increasing in absolute value as ionization decreases. 

Figure 6a and 6b show the values of log Fig. 

6c shows the value of Again, Figure 6a and 6b show an 

almost constant value inversely proportional to the pressure, 
except for the case of higher pressure and smaller radiation 
intensity. Figure 6c displays a nearly zero value when 
temperature increases (in this case it implies both ionization 
and radiation intensity increases), and negative values nearly 
inversely proportional to the pressure values when temperature 
decreases in the grid. 

Figure 7a and 7b exhibit the quotient ( )» arid Figure 

s 

7c shows the ratio (^r/^r )• Note that this coefficient displays 
higher sensitivity to the inelastic collisions since it gives 

« 

greater weight to high velocity atoms. It is expected to be the 
one that reflects the most "noise" in the calculations of the 
antisymmetr ical parts of the atcxn distribution function. 

Figure 7a and 7b show a decreasing ratio with temperature 
and a slope which increases with decreasing pressure. Figure 7c 
shows a negative value at low ionization nearly inversely 
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proportional to pressure, a maximum value, and slowly decreasing 
values at higher temperatures. 

The first-order corrections to the radiative transfer 
equations can be expressed in terms of the functions f ® , since 
schematically 


^ * ''jS' atoms 

and 


= S® r flf. •) + 1 

'■''jS -'electrons '■ ^s ^protons J * 


In these formulas, first-order corrections apply for small 
values of the thermodynamic forces. In the case of high 
radiation flux (for instance, in the winds of hot stars), the 
corrections can rise to zero order, giving very asymmetrical line 
profiles even in the fluid frame. 

Of course, when velocity reaches certain limits depending on 
thermodynamic force values, the method described gives values 
of fj which cannot be considered small first-order quantities. 

At higher velocities and for some angles, our method can give the 
unphysical result f^ > F. In such cases, some cutoff procedures 
are available (Shvarts et al . 1981), which can be introduced in 
out method. However, if calculations are needed for the 
realistic distribution function at high velocities, more 
iterations must be performed on the Newton-Raphson technique, and 


thus the linear relation between the fluxes and the 


35 


thermodynamical forces would not apply. This procedure also 
fails when locally defined coefficients do not apply, reflecting 
the departure fron p-LTE. 


VII. CONCLUSIONS 

For the charged particle collision terms, we used only the 
Landau term in most calculations, but in some of the later 
computations the Boltzmann term was introduced to account for 
collisions between charged particles that result in large angles 
of deflection. Frcm our results, including the Boltzmann term 
leads to negligible corrections in the principal range of 
particle velocities. 

One limitation of the computational procedure used here is 
that relativistic corrections were not accounted for in our 
treatment of particle-photon elastic collisions. Such 
corrections can affect extreme cases. 

From our results, we conclude that it is a fairly good 
approximation to assume f^(u) proportional to p. However, a 

a 

polynomial expansion of fj^(v) is quite impractical, and we obtain 
much better numerical behavior by using a finite partition in v- 

I 

space . 

We also conclude that some of the usual theoretical models 
and the interpretations of observations have to be revised since 
they have underestimated or neglected important transport 
phenomena. One of these is the radiation- induced flow in some 
stellar atmospheres vdiere this flow can play an important role in 
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the energy balance. Another such case is the solar upper 
chromosphere and lower transition region where the appropriate 
values of thermal conductivity can solve some important existing 
puzzles. For example, we show in Figure 9 the temperature 
structure of the region according to model C of Vernazza, Avrett, 
and Loeser (1981). In the region vrtiere the steep temperature 
rise starts, we have plotted the conductive flux corresponding to 
our (p = 0.1 dyne cm“^, W = 0.0022 ,Tj^ = 8000 K) calculated 
conductivity and the one according to Spitzer's (1962) 
coefficient. In Fig. 9c we have plotted the conductive flux 
divergence corresponding to both fluxes. It is striking when 
comparing this figure with Figure 49 from Vernazza, Avrett, and 
Loeser that we get flux divergence in this region of the otder of 
the radiative loses, but when using Spitzer's formula, the flux 
divergence is around an order of magnitude smaller. We believe 
that proper consideration of transport phenomena can probably 
explain guite simply the energy balance in this region where all 
sophisticated energy dissipation explanations have failed (the 
same can also be true for prominence energy balance; see Fon tenia 
and Rovira 1985). 

Another example of the importance of adequate theory for 

« 

transport phenomena is the energy, impulse, and matter transport 
across the magnetic fields through atoms (involving ambipolar 
diffusion), but a consistent theory for this is beyond our paper 
and would require the consideration of the dependence of the 

a 

function fj^ on the angle (i.e,, on the full three dimensions of 
impulse phase-space). 
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FIGURE CAPTIONS 
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FIG. 


FIG. 


FIG. 


la. - Log {a/a ) for solar-type radiation field (W=0.0022, 

s 

Tj^=8000 K) for log p = 0, -1,-2. 

lb. - Log {a/a ) for enhanced solar-type radiation field 

s 

(W=0.022, Tj^=8000 K) and log p = 0, -1, -2. The 
curve labeled s corresponds to W=0. 1, Tj^=8000 K, and 
log p = 1. 

l c. - Log (a/Cg) for LTE-type radiation field (W-1, Tj^=T) 

and log p = 0, 1, 2, 3. 

2. - Same as Fig. 1 except for the ratio (r^/r^ ). 

s 

3. - Same as Fig. 1 except for the ratio (\/X_). 

s 

4. - Same as Fig. 1 except for the ratio (D^/D,p ). 

■^s 

5. - Ratio (r^/r^ ) for solar-type radiation field 

'^s 

(W=0.0022, Tp=8000 K) (part a)/ for enhanced solar- 
type radiation field (W=0. 022 , Tj^=8000 K) (part 
b). Labeled s curve corresponds to p=10 dyne cm , 
W=0. l,Tj^=8000 K. Ratio (rj^/rj^ ) for LTE-typ>e 
radiation field (W=l, Tj^=T) (part c). Curve label 
identifying the logarithm of pressure. 

. - Same as Fig. 5 except for the ratios (D^/D^ ) and 

^ °r/°r_ ^ • 
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FIG. 6 




FIG. 7. 


FIG. 8. 


FIG. 9. 


Same as Fig. 5 except for the ratios ) and 

(Xr/Xr ). 

s 

Quotient (f?/pF2 ), for electrons for cases with p=10 
l a 

dyne cm“^ and T=5000 K, W=l, Tr=T (a); T=1 0,000 K, 
W=l, Tr=T (b); T=50,000 K, W=l, Tr=T (c); and 
T=50,000, W=0.1, Tr=8000 K (d). 

Conductive heat flux in Vernazza, Avrett, and Loeser, 
Model C. Part a shows, for a part of that model, the 
temperature T (in K) as a function of the logarithm 
of mass column (in g cm ^). Part b shows the run of 
the logarithm of the heat flux (in erg cm ^ s -^) for 
the usual SH formula and from the actual calculations 
(with solar-type radiation field and pressure p=0. 1 
dyne cm ) . In part c the logarithm of the flux 
divergence is plotted from the shown conductive 
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TABLE 1 

Parameters Used In Calculations 


I 

Case 

I P 

I 

W 

I 

t 

I 

I 

ISC 

I IE-2 

I 

2.2E-3 

I 

8E3 

I 

I 

2SC 

I IE-1 

I 

2.2E-3 

I 

8E3 

I 

I 

3SC 

1 1E0 

I 

2.2E-3 

1 

8E3 

I 

I 

1SE 

I IE-2 

I 

2.2E-2 

I 

8E3 

I 

I 

2SE 

I IE-1 

I 

2.2E-2 

I 

8E3 

I 

I 

3SE 

I 1E0 

I 

2.2E-2 

I 

8E3 

I 

I 

4FM 

I 1E1 

1 

1.0E-1 

I 

8E3 

I 

I 

3ET 

1 1E0 

1 

1.0E40 

I 

T 

I 

I 

4ET 

I 1E1 

I 

1.0E40 

I 

T 

I 

I 

SET 

I 1E2 

I 

1.0E40 

I 

T 

I 


I 6ET 


1E3 


1.eE-t« I 





















APPENDIX 

TABLES OF COEFFICIENTS 


TABLE 2 

Electrical Conductivity 


T 

ISC 

2SC 

3SC 

1SE 

2SE 

s.eeE+e3 

3.33E+12 

3.56E-H2 

3.71E+12 

3.37E+12 

3.71E+12 

e.eoE-i-es 

4.24E+12 

4.54E+12 

4.72E+12 

4.29E412 

4.72E+12 

8.eeE+e3 

6.24E+12 

6.67E+12 

6.90E4-12 

6.29E+12 

6.90E4-12 

1 .eeE-t^04 

8.42E+12 

9.01E-I-12 

9.28E4-12 

8.47E+12 

9.28E412 

1.50E+04 

1.46E+13 

1.56E+13 

1 .60E-t-13 

1.46E4-13 

1.59E413 

2.00E404 

2.15E+13 

2.30E-t-13 

2.36E+13 

2.16E413 

2.35E+13 

2.S0E404 

2.91E413 

3.12E413 

3.19E4-13 

2.92C413 

3.17E+13 

3.00E404 

3.73E+13 

4.00E4-13 

4.10E+13 

3.74E+13 

4.05E4-13 

4.00E404 

5.53E+13 

S.92E413 

6.10E413 

5.54E+13 

5.99E413 

5.00E+04 

7.52E+13 

8.04E-t-13 

8.32E+13 

7.52E+13 

8.11E-I-13 




TABLE 

3 




Thermoelectric 

Coefficient 


T 

ISC 

2SC 

3SC 

1SE 

2SE 

5.00E-f03 

3.11E406 

4.39E406 

1.97E+07 

2.76E+06 

3.12E406 

6.00E403 

3.64E406 

4.99E406 

1.72E+07 

3.29E406 

3.65E-f06 

8.00E403 

4.70E+06 

6.18E406 

1 .61E407 

4.36E406 

4.71E+06 

1.00E404 

5.75E+06 

7.34E406 

1.63E+07 

S.43E406 

5.77E+06 

1.S0E+04 

8.39E406 

1.01E407 

1.85E+07 

8.11E406 

8.41E406 

2.00E404 

1.10E+07 

1 .28E407 

2.12E+07 

1.08E407 

1.11E+07 

2.S0E404 

1.37E+07 

1.55E+07 

2.41E407 

1.35E+07 

1.37E+07 

3.00E404 

1.64E+07 

1.82E407 

2.69E407 

1.61E407 

1.64E+07 

4.00E404 

2.17E+07 

2.34E+07 

3.26E407 

2.15E+07 

2.17E+07 

S.00E-f04 

2.70E+07 

2.87E+07 

3.82E407 

2.68E407 

2.71E+07 


3SE 

4nirl 

3ET 

4ET 

SET 

6ET 

4.00E4-12 

4.45E+12 

1.96E-t-12 

1.10E412 

4.71E-t-11 

1.75E+11 

S.09E+12 

5.63E+12 

4.51E+12 

4.26E412 

3.29E412 

1.87E-I-12 

7.43E+12 

8.19E4-12 

7.77E+12 

8.66E-I-12 

9.36E412 

9.41E-H2 

9.99E+12 

1.10E+13 

1.04E413 

1.17E+13 

1.34E413 

1.53E-I-13 

1.72E+13 

1.87E+13 

1.77E+13 

1.98E413 

2.24E413 

2.59E-«-13 

2.52E+13 

2.75E+13 

2.59E-f13 

2.68E413 

3.24E413 

3.70E413 

3.41E413 

3.70E4f3 

3.48E413 

3.85E-H3 

4.32E+13 

4.90E413 

4.36E4-13 

4.73E+13 

4.44E+13 

4.90E413 

5.47E+13 

6.18E-H3 

6.4SE-I-13 

6.98E-I-13 

6.53E-H3 

7.18E+13 

7.96E+13 

8.93E-«-13 

8.73E413 

9.45E+13 

8.B3E413 

9.66E413 

1.07E414 

1.19E+14 


3SE 

4FM 

3ET 

4ET 

SET 

6CT 

4.40E406 

S.53E406 

6.67E405 

4.81E405 

-S.83E404 -4.32E405 

5.00E+06 

6.13E406 

-4.52E+07 

-5.74E+05 

1.06E406 

B.79E405 

6.19E406 

7.36E406 

4.44E406 

5.05E406 

7.37E+06 

3.25E+06 

7.35E+06 

8.57E406 

S.42E406 

5.46E406 

S.60E406 

6.44E-f06 

1.01E+07 

1.15E+07 

8.11E406 

8.14E406 

8.17E406 

B.22E-f06 

1.29E+07 

1.44E407 

1.08E+07 

1.08E407 

1.09E+07 

1.09E407 

1.55E+07 

1.71E+07 

1.35E+07 

1.35E+07 

1.36E407 

1.36E407 

1.82E+07 

1 .99E407 

1.62E407 

1.62E407 

1.63E-f07 

1.63E407 

2.3SE407 

2.53E407 

2.15E+07 

2.16E+07 

2.17E+07 

2.17E+07 

2.87E+07 

3.06E+07 

2.69E407 

2.70E+07 

2.70E+07 

2.71E+07 



TABLE 4 

Thermal Conductivity 


T 

ISC 

2SC 

3SC 

1SE 

2SE 

5.008403 

3.308404 

6.168404 

1.148405 

9.598403 

3.338404 

6.008403 

3.468404 

6.728404 

1.138405 

1.058404 

3.518404 

8.008403 

3.868404 

7.808404 

1.218405 

1.388404 

3.948404 

1 .008404 

4.388404 

8.958404 

1.348405 

1.908404 

4.528404 

1 .508404 

6.468404 

1.248405 

1.768405 

4.058404 

6.808404 

2.008404 

9.848404 

1.698405 

2.318405 

7.558404 

1.058405 

2.508404 

1 .478405 

2.298405 

3.028405 

1.258405 

1.578405 

3.008404 

2.118405 

3.048405 

3.898405 

1.918405 

2.278405 

4.008404 

3.938405 

5.068405 

6.208405 

3.748405 

4.238405 

5.008404 

6.508405 

7.878405 

9.348405 

6.328405 

6.998405 


TABLE 5 

Thermoambipolar Coefficient 


T 

ISC 

2SC 

3SC 

1SE 

2SE 

5.008403 

9.448410 

1.468410 

9.558409 

8.938410 

9.448409 

6.008403 

1.288411 

1.848410 

9.208409 

1.248411 

1.288410 

8.008403 

2.108411 

2.728410 

9.948409 

2.098411 

2.108410 

1.008404 

3.138411 

3.778410 

1.138410 

3.158411 

3.138410 

1 .508404 

6,628411 

7.118410 

1.598410 

6.748411 

6.628410 

2.008404 

1.148412 

1.158411 

2.148410 

1,178412 

1.148411 

2.508404 

1.768412 

1.718411 

2.768410 

1.608412 

1.768411 

3.008404 

2.518412 

2.398411 

3.478410 

2.568412 

2.518411 

4.008404 

4.438412 

4.118411 

5.098410 

4.518412 

4.438411 

5.008404 

6.898412 

6.348411 

7.018410 

7.008412 

6.898411 


3SE 4FM 

6.22E404 7.19E4B4 
6.77E+«4 7.84E+B4 
7.SeE+e4 9.19E404 
9.12E494 l.eSE+95 
i.28E4e5 i.seE^es 
1.77E+85 2.07E+85 
2.41E40S 2.81E+05 
3.22E40S 3.73E+05 

5.418485 6.18E+85 
8.44E48S 9.53E485 


3ET 4ET 

7.39E484 7.31E484 
6.84E484 8.14E484 
2.18E484 5.54E404 
1.74E484 2.12E484 
4.41E404 4.92E404 
8.61E404 9.55E404 
1 . 45E405 1 . 60E405 

2.22E405 2.44E405 

4.36E405 4.78E405 
7.37E+05 8.05E+05 


5ET 6ET 

7.24E404 7.21E+04 
8.28E404 8.12E404 
9.35E-I-04 8.72E+04 
3.62E404 7.94E404 
5.56E404 6.43E404 
1.07E405 1.22E40S 

1,798405 2.02E405 

2.728405 3.068405 

5.288405 5.918405 
8.878405 9.878405 


3SE 4FM 

1.468409 2.048408 

1.848409 2.488408 

2.738409 3.478408 

3.778409 4.608408 

7.118409 8.098408 

1.158410 1.258409 

1.718410 1.798409 

2.398410 2.448409 

4.118410 4.058409 
6.348410 6.118409 


3ET 4ET 

-5,758408 -2.838407 
4.528408 -1.868408 
2.088409 2.188408 
3.178409 3.168408 

7.028409 6.898408 

1.208410 1.218409 

1.258410 1.838409 

6.248410 2.448409 

5.948410 8.688409 

6.498410 8.888409 


5ET 6ET 

•1 .348406 -1 .048405 
-7,768406 -3.648405 
3.488407 5.088405 
3.138407 3.298406 

6.798407 6.778406 

1.208408 1.188407 

1.658408 1.838407 

2.618408 2.628407 

4.218408 4.528407 

4.058408 6.788407 



TABLE 6 

Photoeloctric Coefficient 


T ISC 2SC 3SC 

5. eeE-feS -4.41E-fB7 -6.79E4«6 2.61E-fB6 

6. eeE-i«3 -5.eSE4e7 -7.30E+B6 2.63E-f06 
8.eeE-f03 -6.33E-f07 -8.34E406 3.e4E-f06 

1 . eeE+04 -7 . 62E^7 -9.41 E-«^6 3 . S8E406 
1 . 50E+04 -1 . 09E+e8 -1 . 22E+07 4 . 87E+06 

2 . eeE+e4 -1 . 42E4«8 -1 . 52E407 5 . 96E-fe6 

2 . 50E+04 -1 . 76E408 -1 . 84E407 6 . 85E+06 

3.00E+04 -2.09E408 -2.1SE-f07 7.S4E+06 

4.00E-f04 -2.76E-f08 -2.80E+07 8.S2E-f06 
5.00E-t-04 -2.93E-»«8 -3.48E407 9.03E406 


1SE 2SE 3SE 4FM 

-3.54E+08 -4.41 E+07 -6.79E406 -2.93E406 
-4.21E408 -S.0SE-f07 -7.30E406 -3.01E-f06 
-5.55E+08 -6.33E+07 -8.34E406 -3.17E+06 
-6.90E408 -7.62E-I«7 -9.41E406 -3.37E^6 
-1.03E409 -1.09E^8 -1.22E407 -4.08E-f06 
-1.37E+09 -1.42E408 -1.52E+07 -5.02E406 
-1.76E+08 -1.84E+07 -6.13E+08 
-2.09E-f08 -2.16E407 -7.38E406 
-2.77E+08 -2. 81 E+07 -1.01 E+07 
-3.47E+07 -1.30E+07 


TABLE 7 

Photoofflbi polar Coefficient 


T 

ISC 

2SC 

3SC 

5.00E+03 

5.93E+10 

5.48E+09 

1.S4E+09 

6.00E+03 

8.29E+10 

7.49E+09 

1.6SE+09 

8.00E+03 

1.41E+11 

1.24E+10 

2.12E+09 

1.00E+04 

2.14E+11 

1.85E+10 

2.75E+09 

1.50E+04 

4.60E+11 

3.92E+10 

4.78E+09 

2.00E+04 

7,96E+11 

6.77E+10 

7.41E+09 

2.50E+04 

1.22E+12 

1.04E+11 

1 .06E+10 

3.00E+04 

1.73E+12 

1 .49E+11 

1.44E+10 

4.00E+04 

3.01E+12 

2.62E+11 

2.35E+10 

S.00E+04 


4.09E+11 

3.S0E+10 


1SE 

2SE 

3SE 

4FM 

6.53E+10 

5.93E+09 

5.49E+08 

5.92E+B7 

9.02E+10 

8.29E+09 

7.50E+08 

7.90E+07 

1.51E+11 

1.41E+10 

1.24E+09 

1.27E+08 

2.24E+11 

2.14E+10 

1.8SE+09 

1.85E+08 

4.62E+11 

4.6OE+10 

3.92E+09 

3.80E+08 

7.78E+11 

7.96E+10 

6.77E+09 

6.43E+08 


1.22E+11 

1.04E+10 

9.77E+08 


1.73E+11 

1.49E+10 

1.38E+09 


3.02E+11 

2.62E+10 

2.42E+09 



4.09E+1O 

3.75E+09 


TABLE 8 

Photothermal Coefficient 


T 

ISC 

2SC 

3SC 

1SE 

2SE 

3SE 

s.eeE483 

1.16E+07 

2.26E+07 

4.16E+07 

-1.26E+07 

1.16E407 

2.24E407 

6.eeE+93 

1.44E+07 

3.13E+07 

4.73E+07 

-2.44E+07 

1.35E+07 

3.12E407 

8.eeE493 

1.51E+07 

5.30E-t«7 

6.58E-t07 

-4.62E+07 

1.90E407 

5.19E+07 

i.eeE494 

1.98E497 

7.74E+07 

9.08E407 

-1.43E408 

1.62E407 

7.75E+07 

1.50E-te4 

2.S1E+06 

1.57E+08 

1 .73E408 

-7.80E+08 

3.S4E407 

1.64E-t08 

2.eeE494 

-1 .22E-f08 

2.63E+08 

2.84E-f08 

-1.97E+09 

-4.01 E+07 

2.SSE408 

2.SeE+94 

-7.05E+07 

3.73E+08 

4.12E+08 


-1.05E+08 

3.87E408 

3.eeE-«^4 

-'3.09E-»^08 

4.01 E^8 

5.86E+08 


-S.80E-f08 

4.91E408 

4.eeE494 

-1.72E+09 

2.73E+08 

8.98E+08 


-4.56E+08 

S.80E+08 

s.e0E+e4 


1.19E-f09 

1.29E409 



1.09E409 


2.32E407 

3.24E-fB7 

5.46E4B7 

8.14E407 

1.67E4B8 

2.88E4e8 

4.19E488 

8.eSE4e8 

8.73E408 

1.67E489 
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